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Abstract 


There are a finite number of identical aircraft each of 
which contain a number of different types of components 
which fail at different rates. In order for an aircraft to 
be operational, all of its components’ must be operational. 
Each component type has a finite number of spares. Failed 
components are repaired at ae single server facility. 
Simulation is used to study the effect of 6 different repair 
scheduling policies. The repair policies are compared on the 
basis of average number of operating aircraft at the end of 
a mission period of one week. It is found that a repair 
policy which first repairs the component of the type with 
the fewest operating components is the best. In particular, 
it is much better than first-in, first-out, and also may 
well improve upon a policy that serves the longest waiting 
line first. A simple spares stockage policy is developed and 


evaluated when the above scheduling policy is in use. 
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I. INTRODUCTION 


This project was defined by and is under the direction 
of Professor Donald P. Gaver of the Department of Operations 
Research, of the Naval Postgraduate School. The objective is 
to study the means by which a maintenance policy might be 
more effective in support of a client. In this study, the 
client is a detachment of U.S. Navy P-3 aircraft. 

This thesis presents a comparison of six different types 


of repair policies which a server may implement on a queuing 


line of different components. These components are 
subsystems of larger, identical systems (eg. radios and 
aircraft). Each of these components are mission-essential to 


the aircraft and fail independently of each other. There are 
Ns Spares allocated to each component i, i =1,...,I1. The 
repair facility and the spare provisioning constitute the 
maintenance effort of the detachment. The repair facility 
has a single server who must decide on priorities for 
repairing the different types of failed components 
(Figure 1). 

Simulation is used to study six repair policies for the 
Server The £irst priority policy is First-In, First-Out; 
the second policy serves the most numerous component type in 
queue first; the third serves first the largest product of 
the numbers of component types and individual component 
traffic intensities; the fourth's criteria is to serve the 
least frequent fallure component first (smallest failure 
rate); the fifth repairs the most frequent failure component 
first (largest failure rate); and the sixth scans the 
current operating inventory for each component and repairs 
the component with the fewest operating. The six policies 
are compared on the basis of the average number of aircraft 


that have all I mission critical components up after a 


mission of length T. The results of the simulations are 


presented in tabular form in Appendix B. 
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Figure 1.1. System Diagram with Repair Facility 


Maintenance and Material Management data from the Naval 
Aviation Logistics Data Agency (3-M/NALDA) are used for 
estimating failure rates and repair rates. In addition, 
extensive interviews with personnel at the Aviation 
Intermediate Maintenance Department at Moffett Field, 
California; and Naval Air Systems Command at Washington 
D.C., supplemented the data. The interviews indicated a need 
to look at the small detachment (less than six aircraft) 
problem as well as the short detachment problem (less than 
90 days). 

The Aviation Detachment consists of a small number of 
aircraft, each of which contains large numbers’ of 
components, sub-components,sub-sub-components, and so forth. 
Bach of these systems contribute to the availability ofa 
given aircraft's readiness for flight. Each aircraft, in 
turn, contributes to the detachment's ability to carry out 
its mission for some duration, which is typically not a long 


period of time. 


Simplifying assumptions are made in order to reduce the 
size and complexity of the analysis, aid in model 
verification, and provide understandable results. The only 
equipment considered for this study were pieces of 
electronics gear which were designated as mission-critical 
and failure-prone. The length of a mission is measured in 
flight hours and approximates one week's worth of flying. 
Resupply of the detachment occurs at that time; all parts 
are replenished. 

The simulation models were implemented using FORTRAN 77 
programs that are listed in Appendix A. The user inputs 
failure rates of specified parts, total mission time before 
resupply, total components, and total initial spares per 
component. The output displays the average number of “up 
aircraft" at the time of mission completion, as well as 
distributional information. A detailed discussion is 


provided in Chapter III. 


II. NATURE OF THE PROBLEM 


The detachment of small numbers of P-3 Orion aircraft 
by the U.S. Navy to remote sites around the world isa 
Current method of operation by Fleet Commanders. These 
detachments are of relatively short duration, about six to 
eight weeks. They are required to be self-sustained or 
sustainable through air resupply. Small and measured amounts 
of replacement components arrive by P-3 or C-130 aircraft. 
The ability of these detachments to perform their mission is 
heavily dependent on (a) spare parts flown in initially and 
then at specified intervals, and (b) maintenance repair 
policy and capability. 

In the case of the first, (a), considerable effort has 
been expended to address the spare part requirements for P-3 
detachments to remote areas. However, the bulk of these 
consider larger detachments than three aircraft [Ref 1], and 
longer duration than six weeks. The case of provisioning 
smaller, shorter duration groups is most frequently achieved 
by senior maintenance personnel using corporate knowledge 
and experience to derive the requirements. In the case of 
the second (b), little information is available on the study 
of the effects of different repair policies. General 
experience by the author indicates that repair policy is a 
function of the current maintenance administration (policy 
as it relates to which part to fix when). 

The overall goal of the Chief of Naval Operations is to 
achieve at least seventy-two percent of fully mission 
Capable aircraft in a squadron. In whole numbers of 
aircraft, this translates to two of three in a small 
detachment. But this goal belies the fact that for a small 
group of planes to accomplish its mission, clearly all the 


aircraft must be available most of the time. 


The equipment selected for the present analysis was 
avionic gear which had, (a) mean time between failures 
short enough to ensure a reasonable chance of failure during 
a period of detachment, (b) was determined to be mission 
essential by proper authorities, and (c) was repairable or 
replaceable by the detached personnel. Each piece of gear 
was considered to fail independently of the others. The 
overall mission availability of an aircraft is modelled as a 
series system. If all the components are up, then the 


aircraft is up. 


Availability 





Figure 1.2. Minimum Path Representation of N 
Component System 


The failure interarrival times at the queue of all 
components are assumed to be exponentially distributed. The 
service times are also assumed to be exponential. For the 
first series of simulation runs, the failure and service 
rates were the same for all parts. For the latter simulation 
runs, individual failure and repair rates for components are 


used. 


III. MODEL DESCRIPTIONS 


A. FAILURE MODELS 
There are a variety of approaches to the modelling of 

component supply in a multi-component system with spares 
allocation. Generally speaking, as components’ fail and 
Spares are utilized for replacements, spare part inventories 
can be entirely depleted. Any additional failures which 
occur beyond that inventory can result in systems as a whole 
being unavailable as they await maintenance. The rate at 
which parts arrive at the service center is a constant until 
that point is reached because the number of parts in use at 
any given moment,based on the operating aircraft, is 
constant. When the number of operating systems (aircraft) 
starts to drop, so do the corresponding arrival rates. A 
model incorporating this effect is referred to as the 
decreasing arrival rate model. This is a simple situation to 
Simulate, and to obtain results for. However, it is not easy 
to derive analytical results for it. A more analytically 
tractable model is to assume arrivals of a part type to the 
service center form a Poisson process. This will be 
hereafter referred to as the constant arrival rate model. 
This latter approach is rationalized as follows; the 
remaining systems are required to increase their work load 
(flying hours) to compensate for the loss of a system(s), 
ergo more strain on the remaining systems. In addition to 
this, though, a higher (constant) arrival rate would 
represent a conservative approach to calculating spare 
requirements. A higher arrival rate would yield greater 
Spares allocation in provisioning (most provisioning models 
use failure rates as the prime method of specifying spare 


requirements). 


Initially, two simulation models were written to compare 
the two models for numbers of failures: constant, and 
alternatively decreaSing, arrival rates. Ideally, if the 
constant rate simulation results did not depart excessively 
from the actual decreasing rate simulation, then it would be 
desirable to use the constant rate model because it could be 
studied analytically. 

Monte Carlo simulations were written in Fortran 77 and 
utilized a proven random number generator (LLRANDOM II). The 
Simulations modeled aircraft components failing and being 
replaced by available spares; failed components’ were not 
serviced by a repair facility. The simulations used event 
stepping from failure to failure with exponential times 
between failures. To facilitate the analytical analysis, the 
number of aircraft was set at three, the total number of 
types of components per plane at two, and the number of 
Spares for each type of component at one. The failure rate 
of component 1 is 0.02 failures per hour; the failure rate 
of component 2 is 0.0143 failures per hour. Each aircraft 
requires both component types to be in operation so that it, 
in turn, may operate. The generated output was the average 
number of up aircraft at the end of a specified time. The 
constant arrival model was verified through calculation 
(Appendix A). This simulation was then modified slightly to 
cause the failure rate to decrease when whole systems 
dropped off line to create the second model. Both 
Simulations had five hundred replications. The results are 
listed below in Table III-1. 

TABLE III-1l. RESULTS OF MONTE CARLO SIMULATIONS OF 


DECREASING, CONSTANT, AND ANALYTICAL 
FAILURE RATE MODELS 


Expected 
Initial Initial Time in Up 

Model Aircraft Spares Hours Aircraft 
Decreasing 3 1 120 e022 
Constant 3 1 120 Om 22 

Analytical 3 al 120 0.0201 


(Constant ) 


The conclusion drawn from the simulations 1s that the 
constant failure model does not describe well the situation 
of the failure rate for the components decreasing as whole 
systems drop off line. However, when whole systems or 
aircraft drop off line for maintenance, the remaining 
aircraft must fly the same work load. This results in 
subjecting the components on the remaining planes to more 
hours in operation. Therefore, the assumption that demand 
rates for spares decrease as whole systems drop off line is 
not entirely correct either. This situation of increased 
hours of operation on remaining aircraft is not easily 
quantifiable and might be consequential. The disparity 
between the models may be less dramatic for this reason. For 
this reason, the constant model is more commonly used in 
application (e.g., in various METRIC models at the RAND 
Corporation) (Ref. 2]. Reality is probably a middle ground 
between the two. However, the constant arrival rate model 
will be used in this with the caveat that it produces a 


conservative approach to provisioning. 


B. QUEUING MODELS 

The repair scheduling policy that a maintenance facility 
uses directly affects the available inventory of spares. How 
dramatically a given policy may influence systems as a whole 
is of concern to any maintenance supervisor. Six repair 
Simulation models were written in Fortran 77 to study the 
effects of different repair policies on keeping whole units 
or systems operating. Model algorithms are listed in 
Appendix C. Example Fortran code is listed in Appendix D. 

1. Simulation Background 

The baseline Simulation models an aircraft 
detachment located at a remote site. The detachment consists 
of several aircraft, the collection of which are serviced by 
a Single maintenance facility. This facility has a single 
server who can work ona Single part at a time. The server 
draws his parts from a single queue. As soon as he completes 
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a repair, he immediately commences work on the next part if 
a backlog exists. Each aircraft contains several component 
types, each of which has a spare allocation (although the 
allocation may be zero). Each aircraft requires one of each 
component to operate. The detachment is expected to conduct 
operations for a period of time, and to be self-sufficient. 
A resupply occurs at the end of this period and all Spares 
are replenished. The primary measure of effectiveness is 
average number o£ operating aircraft at the time of 
resupply; variance of the number of operating aircraft is 
also tabulated; the empirical distribution could be 
tabulated if desired. For this study, resupply occurs at the 
end of one week or approximately one hundred and twenty 
meLgnt hours. 


The assumptions which were incorporated in the 
Simulation are as follows: the interarrival times of 
component failures are described by an exponential 


distribution; each component haS a unique failure rate and 
fails independently of other components; failed components 
are instantaneously removed from the aircraft, replaced with 
available spares, and placed in either the repair queue or 
directly with the repair server; the arrival of failed 
components to the repair facility is described by a Poisson 
process with a rate equal to the sum of all component 
failure rates times the number of detached aircraft; 
cannibalization or the interchange of components to maximize 
working whole aircraft will occur; an aircraft needs all of 
its components to operate and anything less constitutes a 
down aircraft. Consequently, arrivals of failed components 
to the repair facility form a Poisson process with constant 
rate which does not change if the number of whole systems 
Change. However, if no aircraft are operating, then the 


arrival rate is Q. 


A component which arrives for repair is of type i with 
probability: 

As 7 (Xa) 2 tet ae Ae 
Component i's service time is described by an exponential 
distribution with parameter Ms (mu); no balking occurs at 
the queue; the queue has potential length equal to the total 
number of components in the system; and the repair facility 
is capable of repairing all components. 

2. imulation Structure 

All simulations evolve from a basic "failure-repair" 
model with a First-In, First-Out service policy. These 
repair models simulate continuous time systems by using 
"event stepping" of a simulation clock time from failure to 
repair, failure to failure, repair to failure, repair to 
repair events. 

All simulations are structured with a main program to 
read input files, generate failures, control simulation 
clock and run parameters, and print statistics; a queue 
subroutine for repairs and policy decision criteria; a 
statistical subroutine to assimilate run data and generate 
output statistics; a random number subroutine to call the 
IMSL library and the LLRANDOM II package to generate arrays 
of pseudorandom uniform variates; and, in all but the 
baseline simulation (First-In, First-Out policy), a priority 
subroutine to assign a component priority based on decision 
criteria. 

Input parameters are total mission duration (in-flight 
hours), total number of components per aircraft (parts), 
total number of spares allocated per component, and total 
number of replications of the simulation. The input file 
contains specific failure and repair rates for each 
component. 

The output from a simulation run contains the arrival 
rate (lambda) at the queue; service rate and traffic 


intensity in those cases where service rates for all parts 


150, 


weme Set Sequal; expected numbers of operating aircraft at 
mission termination, variance and standard deviation; a 
tally of the number of Up Aircraft at the end of each run, 
for all runs; average maximum queue size; and average wait, 
in hours, by a part in the queue. 

The repair policies of the various Simulations differ at 
the waiting line of the queue. Policies adjust part 
positions in line for service by current repair decision 
criteria. 


3. Simulation Techniques 
Pseudorandom number generation iS accomplished by 


utilizing a proven pseudorandom number package resident in 
the computer library titled LLRANDOM II; it was developed by 
Dr. P.A.W. Lewis of the Naval Postgraduate School (to check 
for the effect of the starting seed) (Ref. 3]. A subroutine 
calls LLRANDOM II three separate times, with three separate 
seeds, to generate three distinct arrays to store the 
uniform random variates. These uniform random variates are 
used to generate failure interarrival times, repair times, 
and determine which part has failed. 

Failure interarrival times and repair times are 
exponential and are calculated using an inverse transform 
method. By incorporating this method into the simulations, 
any of a variety of distributions could be selected for the 
repair times (using the memoryless property of the 
exponential distribution heavily in generating failure 
times). The algorithm states that given a cumulative 
distribution function, F(x) = P(Xsx). F(X) is uniformly 
distributed over the interval zero to one. By equating a 
uniform random variate to the CDF, F(X), and solving for the 
inverse of F(X), random variates with the distribution of 
F(X) can be generated. [Ref. 4] 


ak 


Sample data on the number of up aircraft at mission 
termination is assimilated in the STATS subroutine. Sample 


statistics are then calculated where: 


Xs = Number of up aircraft at mission's end on run j 
n = Total number of runs 

Qs = Maximum queue size during simulation run i 

W. = Total amount of time (hrs) spent waiting for 


maintenance by the it” arrival at the queue 


m = Total number of parts which waited in queue for all 
runs 
a n 
X = Efnumber of up aircraft} = Ear /on 
j= 
nN = 
Var{X} = Variance of X = eee =SR)- / (ieee) 
1= 
Sic Ble = Standard Deviation of X = (Var X):? 
m 
Wa = Average wait in the queue = (Gyles / m 
1= 
- n 
Q = Average maximum queue size = 2 / on 
L1= 


4. Simulation Walkthrough 


A simulation run begins with the initialization of 
mission parameters such as _ stock levels of components 
(components in use plus spares), mission clock, subroutine 
parameters for the queue, various flags and counters. The 
initial failure interarrival time is calculated and compared 
to an initial repair time (initially set high). The lesser 
of the two times determines the next event and the mission 
clock is advanced to that exact moment. Stock levels are 
adjusted when the events occur; repairs increment, failures 
decrement stocks. In the case of failures, after the 
appropriate stock is reduced, a component is instantly 
placed in the repair facility where it may either be served 
immediately or join a queue (repair policy will determine 


what position in line the part assumes). After a failed part 
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has arrived at the repair facility, all component stock 
levels are scanned to ensure at least one aircraft remains 
in operation. If all aircraft are down, failures cease, the 
mission clock is advanced to the next repair time, and the 
repair event occurs. If at least one aircraft is up, failure 
and repair events continue as before. In the case of a 
repair event, after the appropriate stock is incremented, 
the next part in line is immediately brought into service 
and a new repair time is calculated. The component stocks 
are scanned again to ensure at least one aircraft is up. If 
the event where previously no aircraft were up but because 
the correct type of part was fixed one now exists, a failure 
interarrival time is then calculated and event comparison 
transpires as before. If the all aircraft are still down 
after a repair, the clock again advances to the next repair 
time and that event occurs. 

The mission clock is allowed to be advanced from event 
to event until the scheduled time exceeds the predetermined 
mission duration. When this value is met or exceeded, the 
Simulation run stops, data is gathered on those left waiting 
for service in the queue (waiting time) and the number of 
operating systems is recorded. The run number is compared to 
the total number of replications value input at the start. 
Pt more runs are required, a mission profile is 
reinitialized and the process repeats itself. Otherwise, 
statistics described in previous sections are computed and 


then displayed (at the terminal). 
5, First-In, First-Out Model 


This model is the basic simulation from which all 
others were extended. The First-In, First-Out or FIFO is a 
common form of maintenance scheduling policy. When a part 
arrives in the queue it is7 placed at the end of the line. 
When the server becomes available to work, the part at the 
head of the line or first to arrive, is served next. There 


are never any interruptions once service starts. 
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6. Dynamic 1 Model 
This policy adopts the following criteria: use a 
First-In-First-Out routine until such time as’ there are two 
or more of a given type of part in line, at which time the 
more numerous part type receives first priority for repair. 
For example, there are three types of parts in the queue 
awaiting maintenance, parts A, 8B, and C. There are two part 
A's, three part B's, and one part C. Parts B would move to 
the front of the line and be serviced first as long as they 
remain the most numerous. 
7. Dynamic 2 Model 
This service rule tells the repairman to count the 
number of parts of each type of component and multiply this 
by the individual traffic intensity of the component, then 
pick the largest value. The rule in the form of equations is 
as follows: 
1 = da / Ma = Traffic intensity of component i 
Priority = Max{ 0. x Ni(t)} 


Where: 
Nia(t) = the number of parts of type i in the service 
line at time t 
Ms = the service rate of component i 
As = the failure rate of component i 


8. Dynamic 3 Model 
This scheduling policy determines service priority 
based on the current available stock of operating parts for 
each type component. After the repairman scans’7 the stock 
levels, he reorders the parts in the service line to favor 
the one with the lowest operational inventory. 
9. Failure Rate Priority-Low 
This model assigns priority of service based on the 
failure rate of the type part. Those with the smallest 
failure rate, or least frequently failing, are given first 
service regardless of their arrival time. When the server 


scans the queue for the next part to £1ix, he compares the 
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individual part lambdas to select the next repair. If 
multiple parts are of the same type, a FIFO policy is used. 
The abbreviation FRP-Low is used throughout the text. 
10. Failure Rate Priority-High 

This model is similar to FRP-Low with the following 
modification; priority is assigned to the most frequently 
failing or highest failure rate part. If multiple numbers of 
the same part are being considered, a FIFO policy is used. 
The abbreviation FRP-High is used throughout the text. 


C. DATA 

The primary source for failure and repair rate data in 
this study is the NALDA (Naval Aviation Logistics Data 
Agency) data base. The reason this source is used is because 
it employs extensive error checking algorithms in its 
database. It is considered the best source of clean data 
(error reduced) for P-3's by the experts at Naval Air 
Systems Command and in the fleet. The information it 
supplied was supplemented by personal interviews’ by the 
author with maintenance supervisors at the AIMD, N.A.S. 
Moffett Field, California. 

The interviews cautioned that even though extensive 
error checking algorithms are employed by the data base, 
they were not error free. Precise values for the Mean Time 
Between Failures (MTBF) and the Mean Time To Repair (MTTR) 
would require extensive data analysis. Reasonable 
approximations, however, could be easily extracted. 
Interviews provided an understanding of the sources of the 
data entries and the means to extract reasonable approximate 
values. 

The focus of this thesis is repair policy evaluation. 
The author sought out values for failure and repair rates 
which would test repair decision criteria with a range of 
values and be representative of mission essential avionics. 
The final values are listed in Table IIJI-2. The complexity 
of the P-3 avionics suite makes it necessary to approximate 

> 


This ls done 


subcomponents, which are 


failure and repair rates for components. 


because a component is composed of 


in turn composed of sub-subcomponents. Data information is 
gathered on the subcomponent level but requirements are 
derived on the whole component level. A separate study can 


be done to determine the component failure and 
whole. This 


Failure rates are derived from the 


repair rates 


as a thesis required only approximate values. 
of  £ Fight 
hours between failures (MFPHBF). The repair rates are derived 
AIMD, or 


intermediate maintenance level. It is important to note that 


mean number 


from dedicated maintenance man hours (DMMH) at the 


the precise numerical values are not important as long as 


the relative magnitudes are correct and are reasonably 


representative of actual values. 


TABLE III-2. FAILURE AND REPAIR RATES 


Component,  Lambdas MFHBE Musa DMMH 
1 0.0200 50.0 1.000 aga 
2 0.0143 70.0 0.417 2.4 
3 0.0180 35 ea! 1.000 1.0 
4 Q70026 389 a0 Oss 53 ce 6 
5 070193 51.0 0.666 digo 
6 0.0219 145.8 1733.3 0.8 
7 0.0010 1000.0 0.200 anal) 
8 0.0110 9029 0.400 Za0 
9 0.0300 Srl seo 2.000 C) As) 

10 0.0060 166.7 0.444 Zu 
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IV. SIMULATION RUNS 


A. BACKGROUND 

The six repair policies were exercised through a series 
of different case studies in which input parameters were 
varied identically for all six policies. The Measure of 
Effectiveness (MOE), average number of aircraft in operation 
at the termination of the mission for each repair policy, 
was then compared. Due to the experimental nature of these 
repair policies, a control environment was needed to provide 
a benchmark of performance. Case studies 4 and 5 were an 
attempt to achieve this effect. Their model inputs were set 
equal across all component types (a more detailed discussion 
Follows). Embellishments of the system operating environment 
were conducted in Case 6 and evolved further in Cases 1 
through 3. These refinements subjected the models to ever 
more variable and complex system parameters. 

The input parameters which were altered were: numbers of 
aircraft (planes); spares allocation for each component 
type; and individual component failure and repair rates. The 
total number of component types contained on each aircraft 
was held constant at 9 for all repair policies and cases. 


All simulations had five hundred replications. 
1. Spares Allocation 


Two different schemes for changing spare part 
stockage levels were devised. The first scheme is referred 
to as the "Spare to Plane Ratio" method. Initial spare 


provisions are set equal for all component types (i.e. 
component A's spares are initialized at 12 parts, as are 
component B's, as are C's, etc). As the name of the method 
indicates, the amount which the spare levels are varied is 
in proportion to the number of aircraft in the detachment. 
This ratio changes from 2, to 1.5, to 1, to 0.5, to 0.1666 
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(1/6). For example, if there are six detachment aircraft 
(Planes), all components Spares are varied from 12 (Spare to 
plane ratio of 2), “to 9° (155), to 6 (1), to 3 (05) ace 
(0.166). These provisioning changes constitute the 
differences between Cases 4a-4e, 5a-5e, and 6a-6be. 

The second scheme is referred to as the "K-Standard 
Deviation” method. Predictions of the average number of 
failures a component will experience during a mission period 
iS a measure of the demand for spares. The Poisson process 
is employed to estimate the mean number of spares each 
component would require during a mission period. The 
sensitivity of this estimate of the demand for spares is 
tested by adding to or subtracting from the mean K standard 
deviations OF Spares (from the appropriate Poisson 
distribution). For example, if the mean number of spares 
plus two standard deviations were stocked, the demand for 
parts could be expected to be met 95% of the time (without 
repair). The provisioning calculations proceed as follows: 

Expected Demand for Part i = fa x t x Aa] = EB(X,) 

Standard Deviation of Demand = ( E[X:i] )°-°? 

Spare Policy Component i = E(Xi] + K x ( E(X:i] )°? 

Where: 

a = Number of aircraft 

t = Total mission length 

\a = Failure rate of Component i 

K = Standard Deviation factor 
The result of this policy for components having different 
failure rates is stock levels which are different for each 
component type. The most frequently failing component 
receives the most spares. Variations of the value of K 
determine the difference in Cases la-lh, 2a-2h, and 3a-3h. 

2. Failure and Repair Rates 

Two separate parameter schedules are used in the 
experimentation process. In the first, failure and repair 


rates are set equal for all component types (1.e., Component 
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A's failure rate is 0.02 failures per hour and its' repair 
rate is 0.84 units per hour, as are Component B's, etc). 
This input schedule is used for Cases 4a-4e, and 5a-5e. The 


second schedule uses individual failure and repair rates 


listed in Table III-2 (Failure rates = Lambda:, Repair rates 
= Mu:). This parameter schedule is used for Cases la-lh, 2a- 
2h, 3a-3h, and 6a-6h. 

3, Output 


The output of each case is given in Appendix B. The 
listed results include the Mean number of operating aircraft 
at mission termination (Average Up Planes), the Variance of 
the number of up aircraft, the Standard Deviation of the 
number of up aircraft, the Standard Error of the Mean, and 
the fraction of Up Planes operating at mission termination. 

The number of operating aircraft at the end of a given 
run is determined by analyzing the number of operating parts 
of component type i, i= 1,...,I. The component with the 
minimum number of operating parts determines the maximum 
number of operating planes. The Mean number of operating 
aircraft is calculated from the numerical average of Up 
Aircraft at the end of all replications (runs). This is the 


common sample mean: 


500 
Mean Up Aircraft = ©£ X: / 500 = X& 


Where: 


X, = Number of Up Aircraft on run i 


The Variance of the run sample and the Standard Deviation 


are: 
500 
Variance = gle - X)7 / 499 = Var(X] 
l= 
Standard Deviation = ( Var[(X] )°3 


The Standard Error of the Mean is given by the relation: 
Standard Error = S_ = ( Var{X] / 500 )°° 


The Fraction of Up Planes operating at mission termination 
is used as a measure of effectiveness for the policies. [t 
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is calculated after each subcase is completed and is derived 
by dividing the Mean number of operating Aircraft (E{X]) by 
the total number of detachment aircraft. (Ref. 5] 
4. Random Seeds and Stability 

In order to establish that the simulations had 
reached stability in sampling variance, five of the six 
models were simulated by varying only the random number 
seeds used in generating failures and repairs*. Four sets of 
seeds were utilized. Two-way Confidence Intervals for the 
Means (Expected Number of Up Aircraft) were then computed at 
a level of 95%. The distribution of the Sample Means was 
assumed to be Normal in accordance with the Central Limit 
Theorem. The results are listed in Appendix B. and 
demonstrate that the sampling variance had stabilized. For 
any given policy, the Means of any combination of three 
random number sets fell within the confidence interval of 
the remaining fourth set. For any given set of random number 
seeds, the six repair policy Means established a_ rank order 
by magnitude. This relative ranking between policies did not 
change for all four seed sets. As a result, it was 
determined that reliable output from the simulations could 


be anticipated with five hundred replications. 


B. CASE DESCRIPTIONS 

1. ‘Ga 4- 

These cases were designed to give a baseline measure of 
how the models would compare in a environment’ where 
interrelational effects between parameters such as failure 
and service rates, and stock levels, could be minimized. 
Case 5 differs from Case 4 in the number of total aircraft. 
Subcases of Case 4 and 5 delineate variations in initial 
spare provisioning as described by the "Spare to Plane 


Ratio" method. Inherent to the background of both cases are 


*The sixth model, Dynamic 3. was added after this 
exercise. Stability of variance for it was assumed based of 
the other models' results. 
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failure and repair rates, and total numbers of components 
per plane (9). Failure and repair rates are fixed at equal 
values (0.02 and 0.84 respectively) for all types of 
components. These rates represent average values’ for the 
listing in Table I[III-2. 
2. Cases 1-3 

In these cases, the failure and service rate for the 
components differ. These cases were designed to analyze how 
sensitive the various models were to a dynamic environment 
where several of the parameters are varied (and may have 
potential interrelated effects). For example, high failure 
rates and high service rates in conjunction with large 
quantities of spares may result in the server dedicating a 
Significant portion of his effort towards ae single type 
part. As before, Cases 1, 2, and 3 differ in numbers of 
detached aircraft. Case l contains three; Case 2 contains 


S1x; Case 3 contains twelve. Subcases delineated spare 


provisioning as described by the "K-Standard Deviation" 
method. Cases 1, 2, and 3 differ in their selection of 
values of K. The range of K is tl to -4. Negative values 


of K are selected to explore the effects of decreasing spare 
levels. The actual spare allocations as a function of case 
number are listed in Table B-3 of Appendix B. An excerpt 


from that Table is given below as an example. 


TABLE IV-1 MODEL PARAMETERS: SPARE ALLOCATION OF 
COMPONENTS FOR SIMULATION CASES 


Case la Case 1b Case lc Case ld 
K-Factor +l 0.0 -0.5 -1.0 
PART 1 10 7 6 5 
PART 2 7 5 4 3 
PART 3 9 6 5 4 
PART 4 2 L 0 0 
PART 5 10 4 6 4 
PART 6 ll 8 6 4 
PART 7 1 0 0 0 
PART 8 6 4 3 2 
PART 9 14 ll 9 8 


Zh 


Unique component failure and service rate are input from 
Table III-2. Total components per aircraft are nine. 

3. Case 6 

This was an intermediate step from Cases 4-5 and Cases 
1-3 for parameter selection. Spare stockage levels employed 
the "Spare to Plane Ratio" method as per Cases 4 and 5. This 
determined the variation in the subcases (6a-6e). Unique 
failure and repair rates are initialized at the values 
listed in Table III-2, and as per Cases 1-3. The number of 


component types per aircraft 1s unchanged at nine. 


ZZ 


V. RESULTS OF CASE STUDIES 


The purpose of the case studies is to compare different 
repair policies using as the primary MOE average number of 
aircraft up at the end of amission of finite length. A 
repair policy substantially more responsive to shortages 
then a simple First-In, First-Out routine should stand out. 
The Dynamic 3 policy consistently achieved this result. 
Conversely, the Failure rate priorities are consistently 
poor performers (relative to the FIFO benchmark). The other 
two Dynamic policies are sensitive to the spare provisioning 
policy. Recall the fact that both Dynamic 1 and Dynamic 2 
set priority for service based on total numbers of component 
types awaiting maintenance. High-failure rate parts, which 
are also provisioned more heavily, will tend to consume the 
available server's capacity. Those parts which are stocked 
at lower levels (based on lower failure rates) are relegated 
to the "back of the line". It is the effect of all the lower 
rate failure items being placed at the end of the line which 
causes the Dynamic 1 and 2 to drop in performance (when 
there are unequal numbers’ of spares for each component). 
There is a cross-over point for these two models relative to 
the FIFO policy, that point being when the spare provisions 
are depleted to the point where stock levels’ begin to 
equalize (K factor s -2.25). More detailed case by case 
discussions are provided in this chapter. 

Case study results are listed in Table B-l, Appendix B. 
The graphical summaries contained in this chapter provide a 
visual means to summarize the effectiveness of each of the 
repair policies relative to each other. Since Cases 4-5 
provide a baseline measure for the models, they will be 


discussed first, followed by Case 6, and then Cases 1-3. 
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A. CASES 4-9: sBOQUALBEALEURE AND Se Aine nse 

The equality of failure rates for each component and 
repair rates for —€a€n @e@omponentw should demonstratemere 
relative capability of @€ach repair policy to replenism@ Siem. 
stock levels and keep airplanes operating. A significant 
difference between policies can be determined by comparing 
the Mean and the Standard Error for each data point from 
Appendix B and calculating Normal confidence intervals for 
the mean (also by recalling that for Normal distributions, 
95% of a population distribution should lie roughly within 2 
standard deviations of the mean). By analyzing the data and 
Graphs (Figures 5-1, 5-2), for Cases 4-5, the Dynamic 
policies (1, 2, 3) are consistently more able to keep 


aircraft operating (stocks replenished). 
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Figure Sil. Fraction of © pee Planes as a Function 
of Spares to Planes Ratio (Case 4). 
The difference between them and the FIFO routine is 


Significant at the point where the spare to plane ratio is 
equal to 1.5. Since failure and service rates are equal for 


ZA 


ieee COMPOMeENES, "Dynamic and Dynamic 2 are in essence 
following the same policy as Dynamic 3; that being, repair 
tne lowest stock level Eirst (or conversely the most 


numerous in the queue). 
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Figure 5.2. Fraction of Operating Planes as a Function 
of Spares to Planes Ratio (Case 


The actual differences in the policies are attributed to 
variation within the simulation. The means for the Dynamic 
policies were within 1.96 standard deviations of each other 
(95 % Confidence Interval). The Failure rate priority 
policies were not considered in Case 5 due to their poor 
performance in Case 4. In Case 4, they demonstrated 
performance which was as much as twenty five percent less 
effective than FIFO. The difference between the Failure rate 
Priority routine and FIFO is because FRP-L and FRP-H assign 
priorities one through nine to each of the parts. Since all 
the rates were identical, the priority assigned was one for 


component one, two for component two, etc. Also, since the 
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rates were identical, there is no difference in component 
priorities in FRP-L and FRP-H. 
B. CASE 6: DIFFERENT FAILURE AND REPAIR RATES; IDENTICAL 
SPARES STOCKING 
These simulations address the issue of how the policies 
respond to variable failure and service rates with a 
constant part stockage. Individual component failure and 
repair rates are given by Table III-2. Again, since the 
numbers of spares is the same for all parts, Dynamic and 
Dynamic 2 still reflect the Dynamic 3 policy (numbers in the 
queue relate to the stock levels which relates to Up 
Planes). The Dynamic policies are still consistently better 


than either FIFO or the Failure Rate policies (Figure 5-3). 
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Figure 5.3. Fraction of O cea d Planes as_a Function 
of Spares to Planes Ratio (Case 6). 


Again, the point of significant departure for the Dynamic 

policies and FIFO occurred at a plane to spare ratio of 1.5. 

The Dynamic policies' means fell outside of a 95% confidence 

interval for the FIFO policy. The differences between the 
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Dynamic policies is attributed to variation in the Mean 
(within a 95% confidence interval). The Failure rate 
policies showed improvement. These two policies are now 
dissimilar because failure and repair rates are different. 
Between the two, FRP-H performs significantly better than 


FRP-L (based on a 95% confidence interval about the means). 


C. CASE 1-3: DIFFBRENT RATES, DIFFERENT STOCKS 

In these cases, we see a Change occur in the abilities 
of Dynamic 1 and Dynamic 2 to replenish stock levels. 
Initial stock levels are no longer equal for all parts. They 
vary for each component based on failure rates. High failure 
rate components are supplied the most spares; low failure 
rate components the least. Since Dynamic 1 and Dynamic 2 
repair the component having the most numerous failed parts 
first, the issue of replacing the most needed operating 
stock is skirted. High-failure items can fill the queue and 
still maintain sufficient operating units on all aircraft. 
Other units may fail; however, if they fail in insufficient 


numbers, they may never be served. For example, in a case in 


which the expected demand (failures) is stocked for 
Component 1 and 8 (K = O for 6 planes), their respective 
Spares are 14 and 8. Component 1's failure rate is 


approximately twice eight's (0.02 vice 0.011). Therefore, 
Component 1's arrive at the repair facility twice as often. 
If there are 14 Component 1S waiting for repair and 13 
Component 8s waiting for repair, then the Dynamic 1 and 2 
policies would choose Component 1 to repair next even though 
there are still 6 Component 1s operating and only 1 
Component 8 operating. This policy would result in only l 
plane operating. A comparison of Case 6 with Case 2 will 
demonstrate what happens when the Spare provisioning policy 
is changed from equal spares for all components (Spare to 
Plane Ratio method) to a failure-based system (K-Standard 
Deviation method). The Dynamic 3 policy is concerned only 
with those components which are in most demand (lowest 
725) 


operating stock level), regardless of the number awaiting 
service. It doesn't waste time on repairing those failed 
parts that have sufficient spare stock levels atyeagaue en 
moment. Dynamic 3 was consistently the best of all policies 
(Figure 5-45-55 ace 


ules) 
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Figure 5.4. Fraction of Operating Planes as a Pllmemimen 
of K Standard Deviations of Parts (Case l). 


Once again, the Failure policies do not to perform well. In 
Case 1 with 3 planes, there is a tight bunching of all the 
policies. This is due to the fact that the server is not 
very busy with 3 planes and as a result, the form of the 
repair policy is not critical. As the number of planes 
increase, the policies diverge; the server becomes less able 


to rejuvenate stocks in the face of higher demands. 
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VI. RESULTS AND CONCLUSIONS 


This thesis has presented a comparison of six repair 
scheduling policies. The Dynamic 3 repair model, which 
schedules service in favor of the component with the lowest 
operating stock, yields a Significant improvement of 
available aircraft at the end of a resupply period when 
compared to a esimple First-In, First-Out policy. A key 
factor in repair policy scheduling is the current inventory 
of available operating components. The other models which 
were studied in this thesis failed to address this issue 
directly. As a result, when spares were stocked in unequal 
numbers (the most likely real world situation), the other 
policies (Dynamic, Dynamic 2, FIFO, FRP-L, FRP-H) did not 
perform as well as Dynamic 3. 

The overall goal of the Chief of Naval Operations is to 
achieve at least seventy-two percent of fully mission 
Capable aircraft. The maintenance support of a squadron and 
the supporting aixrwing are responsible for accomplishing 
that goal. It is clear from the studies conducted in this 
thesis that significant differences in aircraft availability 
can result simply by the manner in which a maintenance 
effort schedules its repairs. A policy such as Dynamic 3 can 
assist in achieving the CNO goal better than FIFO. Spare 
stock levels are maintained more effectively regardless of 
the provisioning policy or the workload in the repair 
facility. However, a simple and sensible stocking policy 
increases effectiveness even more; the K-policy illustrates 
this fact. Future study into other repair and stockage 
policies is recommended to study their effects on system 
availability. Other repair policies may look into the 
effects of different distributional assumptions for failure 


and repair times, as well as repair scheduling. It is also 
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necessary to consider the effect of breakdowns of the repair 
facility itself: How can these be best accommodated? 
Finally, it is desirable to reduce the computational effort 
needed to evaluate repair and stockage policies either by 
improving simulation efficiency by Monte Carlo "swindle", or 
by their replacement by analytic approximations. These steps 


are under current examination. 


SL 


APPENDIX A 


A. THE PROBABILITY CALCULATIONS OF UP AIRCRAFT 
The following calculations are used to verify the 


constant arrival rate simulation model with no repair 


FaCillrey.: 
Parameters: 
Ai;(t) = The number of type i components available 
and up during the interval (0,t] 
Fi(t) = The number of type i components’7 which fail 
during the interval (0,t] 
I = The total number of types of components 
a = The total number of aircraft 
ns = The total number of Spares of type i component 
As = The failure rate for component i 


Fach aircraft requires one of each type of component to 
operate. The maximum number of up aircraft at any given 
moment t iS equal to the lowest level of up components for 
any given type of part. Therefore, the number of up aircraft 


at any given time t is: 


A(t) = Min{Aa(t), A2(t), .. ., Axr(t),...,a} (A. 
If we use the constant failure rate assumption, 
component i will fail at constant rate equal toa a, even 


if the number of up aircraft is less than a. We assume an 
exponential failure distribution and cannibalization of 
parts. This leads to approximating the failure time of part 
l1 by a Poisson arrival process. We are interested in the 
probability of a minimum number of aircraft being up at 
time t: 

P{Min # A/C up 2 k}..... 


P{A(t)2 k} 


I 
m P{ As(t)2 k} (A. 2) 
i=1 


P{ # A/C avail. > k} 
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Since different components are assumed to fail independently 
and based on the number of aircraft hours flown (0,t) 

The probability that the number of components that fail 
of type iin the interval (0,t] is less than or equal to 
some fixed value j is: 

P{# parts that fail type i in (0,t]) < j} = 

e-aAie (ad at)* 
ak 


Cu. 


P{F.(t)<j} = , 1f£ j < atna (A.3) 


K 


Stes t)siy = it 7 eee) > a trie 


Now, recall that the number of failed components plus 


the number of available components equals the _ total 
components. 
Fi(t) + Aa(t) = a + ma @ Aa(t) = a + ma —- Fa(t) (A. 4) 


We would like to solve for the probability that the number 
of available components of type 1 at time t is greater than 
Or equal to some number k or P{Ai(t) 2 k}. So substituting 


A.4. for As(t), we solve: 


A 


~Pfa + na - Fa(t) 2 k} ~@ P{Pa(t) < a + ma - kK} = 
P{A,s(t) 2 k}. 


If we let the number j = a+n:a - k and plug into our 
previous solution (A.3) for the number of components that 
fail of type i in the interval (0,t): P{Fia(t) < j} 
atna-k e-aAjit (aAat)® 
P{PF.i(t) < a+t+na- k} = & SuEEEEIeEnn aan (A.5) 
L=0 L! 
= P{A:(t) 2 kK} 
We now solve for the boundary conditions: 
lower: j < a tna, J = at na - kK wk 21 
upper: k < a , since you can not exceed the total number 
of aircraft in the system as an upper Dound. 
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In summary, the probability that the number of up 


aircraft is greater than some fixed number k is given by... 


I atnsa-k e-AAre (a Ast)® 
P{A(t) 2 kK} = E —. if l<ksa 
i=l L=0 L! 
P{A(t) 2 O} = 1 ... the entire probability space 
P{A(t) 2 atl} = 0 ... you cannot exceed the total number of 


aircraft. 

The expected number of up and operating systems at the 
end of some period t is given by the product of the 
probability of a given value and that value: 

E{A(t)} = OxP{A(t)=O} + ILxP{A(t)=1} + 2xP{A(t)=2} + 

The variance of the expected number of up and operating 
systems at the end of period t is given by: 

B{A(t)?7} = OxP{A(t)=O} + 1?7xP{A(t)=1} + 27xP{A(t)} + 
Var{A(t)} = E{A(t)7} - E{A(t)}? 

where: P{A(t)=0} = P{A(t)20} - P{A(t)21} 

The constant arrival simulation without repair is 
verified with the parameters: 

I = 2, a = 3, ma = 1, t = 120, 2 = 1/50, 2 = 1/70 
The probabilities then follow: 

P{A(t)=1} = P{Aa(t)=1} x P{A2(t)21} 


+ P{Aa(t)21} x P{Az2(t)=1} = 0.01479 
P{A(t)=2} = P{Aa(t)=2} x P{A2(t)22} 

+ P{Aa(t)22} x P{A2(t)=2} = 0.00273 
P{A(t)=3} = P{Aa(t)=3} x P{A2(t)23} 

+ P{Aa(t)23} x P{Az2(t)=3} = 0.00001 
E{A(t)} = 1x0.01479 + 2x0.00273 +3x0.00001 
E{A(t)} = 0.0201 , E{A(t)}2 = 0.0004 


E{A(t)?} = 0.0258 
Var{A(t)} = 0.0254 
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APPENDIX B 


TABLE B-l MODEL STATISTICS: EXPECTED NUMBER OF UP 
AIRCRAFT VS SIMULATION CASE NUMBER 
Case la Case lb Case lc Case ld 

Planes & 3 3 3 
K-Factor2 1 0 -0.5 —-j 
Parts/Plane 9 9 9 9 
Lambda/Mu ----Individual Failure, Service rates used----- 
Dynamic lL 

Mean 3.00 2.99 2-93 2.92 
Variance 0.00 0.01 0.07 ORO) ¢ 
S.Dev. 0.00 0.10 Ore27 O27 
Sa 0.00 0.00 0.01 0.01 
%* Up A/C 1.00 0.99 0.98 0.97 
pynamic 2 3.00 2.98 2.89 2.89 
Variance 0.00 0.02 0.09 0.09 
S.Dev. 0.00 0.14 Os32 0.32 
S_ 0.00 0.01 0.01 0.01 
% Up A/C 1.00 0.99 0.96 0.96 
Dynamic 3 
Mean 3.00 Paes) 4) Zao 3 2393 
Variance 0.00 0.03 0.06 0.07 
S.Dev. 0.00 On 17 0.25 0.26 
S_ 0.00 0.01 0.01 0.0l 
% Up A/C 1.00 0.99 0.98 0.98 
FIFO 
Mean 3.00 2.97 2.93 2.92 
Variance 0.00 0.03 0.07 0.07 
S.Dev. 0.00 Oey 0.27 On 7 
S_ 0.00 0.01 0.01 0.01 
% Up A/C 1.00 0.99 0.98 0.97 
FRP-Low 

ean 3.00 2.96 2.91 2.91 
Variance 0.00 0.04 0.08 0.08 
S.Dev. 0.00 0.19 0.29 0.29 
ee 0.00 0.01 0.01 0.01 
% Up A/C 1.00 0.99 0.97 0.97 
FRP-High 

ean 2.99 2.96 2.92 2.92 
Variance oO. 01 0.04 0.09 0.08 
S.Dev. 0.10 0.20 0.29 0.29 
S_ 0.00 0.01 0.01 0.01 
% Up A/C 0.99 0.99 0.97 0.97 


7K-Standard Deviation method used, see Chapter IV 
37S§_ = {(Variance/500)}-*7 = Standard Error 
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TABLE B-1 MODEL STATISTICS: EXPECTED NUMBER OF UP 
AIRCRAFT VS SIMULATION CASE NUMBER* 


Case le Case lf Case lg Case lh 
Planes 3 3 3 3 
K-Factor el =F) 0 -2.5 —3au 
Parts/Plane 9 oe) 2, 9 
Lambda = “tame. -- = Individual Failure rates used-------- 
MO ss ee Individual Service rates used-——_ aaa 


Dynamic 1] 
Mean 


2.89 2.65 2.03 1.93 

Variance 0.12 0.32 0.38 0.39 
S.Dev. 0.34 0.57 0.61 0.63 
s 0.02 0.03 0.03 0.03 
%* Up A/C 0.96 0.88 0.68 0.64 
Dynamic 2 

Mean 2.87 2.63 2.09 1.89 
Variance 0.12 0.33 0.39 0.31 
S.Dev. 0.34 0.56 0.62 0.63 
S_ 0.03 0.03 0.03 0.02 
%* Up A/C 0.96 0.88 0.69 0.63 
Dynamic 3 

Mean 2.91 2.68 2.04 1.91 
Variance 0.09 0.29 0.42 0.36 
S.Dev. 0.31 0.55 0.64 0.59 
S_ 0.01 0.02 0.03 0.02 
%* Up A/C 0.97 0.89 0.68 0.63 
FIFO 
Mean 2.89 2.64 1.96 1.92 
Variance 0.11 0.34 0.38 0.36 
S.Dev. Ones 0.59 0.62 0.60 
a 0.02 0.03 0.03 0203 
%* Up A/C 0.96 0.88 0.65 0.63 
FRP-Low 

ean 2.86 2205 2.02 1.93 
Variance Osi] O27 0.39 0.3333 
S.Dev. 0.41 0.61 0.63 0.57 
S_ 0.02 0.03 0.03 Ons 
%* Up A/C 0.95 0.88 0.67 0.64 
FRP-High 
Mean 2.86 2.65 2705 1.93 
Variance 0.16 0.35 0.38 0.34 
S.Dev. 0.39 0.59 0.62 0.59 
So 0.02 0.03 0.03 0.03 
%* Up A/C 0.95 0.88 0.68 0.64 
“Note: None of the policies differ much, either 


practically or statistically. 
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TABLE B-1l MODEL STATISTICS: EXPECTED NUMBER OF UP 


Planes 
K-Factor 
Parts/Plane 
Lambda 

Mu 


Dynamic 1 
Mean 
Varlance 
oe 


<<, 


%* Up A/C 


Dynamic 2 
Mean 
Variance 
7 


——> 


%* Up A/C 


Dynamic 3° 
Mean 
Variance 
S.Dev. 


% Up A/C 
FIFO 
Mean 


Variance 
7. oY: 


% Up A/C 


a 


Variance 
apc. 


% Up A/C 
FRP-High 
Mean 

Variance 


S.Dev. 
S 


% Up A/C 


=Note: 


AIRCRAFT VS SIMULATION CASE NUMBER 


Case 2a Case 2b Case 2¢ Case 2d 
6 6 6 6 
1. 0 -1 = 
g 3 9 


9 
pls a oo eels Tidiatclduguianlure rates used>——--=—-- 
a hiceivratial servrce rates used-------- 


3599 D935 Dege | Bie. 2 3 
0.01 0.09 Oi 2 O72 
O20 1 2 0.84 0.85 
0.00 0.01 0.04 0.04 
0.99 0.93 0.88 0.87 
5.96 DoS 0:3 a. 0'3 
0.01 O25 i093 09 
0.11 0.49 1.05 iba Ibe, 
0.00 0.02 0.05 0.05 
0.99 0297 0.84 0.84 
6.00 6.00 5.88 oes 8, 
0.00 0.00 QO.11 0.19 
0.00 0.00 0.34 0.43 
0.00 0.00 0.01 0.02 
1.00 1.00 0.98 0.97 
ieee abs) S, a OL D2 3 
0.01 0.02 0.42 0.65 
0.10 0.15 0.65 0.80 
0.00 0.01 0.03 0.04 
0.99 0.99 Webs 0.92 
2394 Disco 4.97 4.58 
0.22 ee 7S), 4.09 
0.46 1.05 a 2 (NOs 
0.02 0.05 0.08 0.09 
0.99 0.95 0.83 0.76 
5.96 abe (fal 4.78 4.69 
0.07 0.66 2.03 2.94 
0.26 0.81 1.42 1 Be, 
0.01 0.04 0.06 0.07 
0.99 0.95 0.80 0.78 


Dynamic 3 looks very good in Cases 2c and 2d. 
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TABLE B-1 MODEL STATISTICS: EXPECTED NUMBER OF UP 
AIRCRAFT VS SIMULATION CASE NUMBER 


Case 5a Case 5b Case 5c Case 5d Case 5e 


Planes pe? 12 12 12 12 
Spares’ 24 18 12 6 2 
Parts/Plane 9 5 3 9 9 
Lambda  9© ses. — Failure rate all components = 0.02----- 
MQ ee Service rate all components = 0.84----- 
jpgehe 

ean 12-300 J a: 2 8.88 S15 0.82 
Variance 0.00 0.16 3.86 S70 2 0.89 
S.Dev. 0.00 0439 P97 1 SPB. 0.3e 
Sa 0.00 0.01 0.09 0.09 0.04 
% Up A/C 700 0.99 0.74 0.26 0.07 
Dynamic 2 
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Variance 0.00 OG. 09 4.07 Sr ae, 0.328 
S.Dev. 0.00 0:29 Ze0Z 194 0.96 
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Dynamic 3 
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% Up A/C P00 O2399 0.73 Ou) Os12 
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Mean oS 10.51 5.61 ie. 25 0.66 
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S.Dev. 07250 2.14 20 3 1.44 0.91 
Sa 0.02 0.10 O13 0.06 0.04 
%* Up A/C 0.99 0.88 0.47 0.10 0.06 


7Spare to Plane Ratio method used for computing spares 
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SPARE ALLOCATION OF COMPONENTS 
K-STANDARD DEVIATION METHOD 
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TABLE B-3 MODEL PARAMETERS: SPARE ALLOCATION OF 
COMPONENTS VS SIMULATION CASE NUMBER 


Case 3a Case 3b Case 3c Case 3d 
+1.0 0.0 -1.0 -2.0 
34 29 23 18 

2 25 21 16 2 

3 31 26 21 16 

4 6 4 Z 0 

5 33 28 23 17 

6 37 31 26 20 

7 3 1 0 0 

8 20 16 12 8 

9 50 43 37 30 
Case 3e Case 3f Case 3q Case 3h 
-2.5 -3.0 -3.5 -4.0 
15 13 10 7 

2 9 7 5 2 

3 13 11 8 6 

4 0 0 0 0 

5 rs 12 9 7 

6 17 15 12 9 

1 0 0 0 0 

8 6 4 2 0 

9 27 23 20 17 
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APPENDIX C 


TABLE C-1. FIFO ALGORITHM 


Key Variable Description; 


Planes: 


MOTSPR : 
TOTMSN : 


mOLP RT: 
TOTRUN : 
MSNTIM: 


LOWLVL : 


STOCK : 


FAILT: 


COMPLT: 


SERVTM: 


QPART: 


QSIZE: 


QTIME: 


The total number of whole systems operational, in 
this case aircraft. 


The initial number of Spares per type part. 

The total mission length in flight hours (120 Hrs.). 
The total number of types of parts. 

The total number of replications or runs (100). 


The mission time running clock, ranges vary from 
0 to 120 hours (TOTMSN). 


An Array which stores the lowest part stock level in 
the event of a given part failure. 


An Array which stores the number of up and working 
components of a given part type (spares + planes). 


A variable which holds the failure time generated by 
a uniform random variable and the inverse transform 
method for an exponential distribution. This is 
compared to the repair time to determine the 
scheduling of next event. 


A variable which holds the repair time, is generated 
like the failure time, and when compared to the 
latter determines the next scheduled event. 


A variable which is assigned a exponential value 
generated by the inverse transform method and is 
added to the current clock to determine the 
completion time of the next repair 


An array which stores component types as’ the wait in 
he queue. 


An array which stores current position a component 
holds in the queue. 


An array which stores the time a part enters the 
queue. 


RHO: The traffic intensity equal to arrival rate divided by 
the service rate in the queue 


SUMLAM: 
SRATE: 


REPAIR: 


FLAG2: 
RUN: A 
Input: 


The sum of all the part failure rates. 

the service rate (Hu) in the queue. 

A counter to track when repairs occur. 

A counter to track if stocks are depleted. 
counter to track run number (0-100) 

Failure rate parameters (Lambda(i)) for each 
component. 
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Repair rate parameters (SRATE(i)) for each component. 


Output: Mean number of up aircraft at missions end 
Variance and standard deviations of above means 
Distribution of up aircraft at missions termination 
(the number of runs with zero aircraft up, the 
number with two, and so forth). 
The average waiting time in the queue 
The average queue line 


Step 1: Initialize simulation parameters: TOTRUN, TOTMSN, 
TOTPRT, SRATE, etc. 


Step 2: Read in repair rates SRATE(i) and failure rates 
LAMBDA(i ) 


Step 3: Call RANDOM; Generate an array to store uniform 
random numbers, store in A( ),AA( ),AAA( ). 


Step 4: Initialize a new mission run: MSNTIM, PLANES, 
SPARE(I), COMPLT, FLAG2, 


STOCK = SPARES + PLANES 
SUMLAM = LAMBDA(1) + LAMBDA(2) + ... + LAMBDA(N) 
RHO = SUMLAM * PLANES / SRATE?? 
Step 5: Verify at least one aircraft is up (FLAG2 = 1) 
Step 6: Generate a failed part: 
FAILT = - ln (U) / (SUMLAM x PLANES) 
This uses the inverse transform method with a 
uniform random number (U) from an IMSL routine 


Se It to yield an exponential interarrival 
ime. 


Step 7: Scan stock levels (STOCK(i)) and verify at least one 
aircraft is up (FLAG2 = 1): 


If STOCK(i) < 1 (FLAG2 = a then check to see if a 
repair has occurred which w 11 create an up 
aircraft. If all aircraft are down, advance to next 
repair time or terminate run. 


Step 9: Test to see if a failure or repair is the next 
event to occur: 


Is FAILT < COMPLT ? I£ so, MSNTIM = FAILT, go to 
step 11; else MSNTIM = COMPLT,go to step 10: Test 
both cases to see if mission has terminated: 

Te MSNTIM greater than TOTMSN ? If so, go to step 


Step 10: Call ALTQ;, Enter queue, compre ce repair and. 
increment appropriate stock. Commence service on 


23, Not calculated when using multiple service rates. 
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Step 


Step 
STEP 


Step 


Step 


Step 


Step 


Step 


Step 


Step 


Step 


oT ie 


12: 


Es: 


14: 


15.: 


16: 


5 ley 


18: 


19: 


20: 


Zi 


the next part if compenents remain waiting 
service, go to step /?. 


Determine which part failed; calculate the 
individual probability a given part i failed: 


Probability Part i failed = LAMBDA(i)/SUMLAM 


Then, generate a discrete CDF with the 

Reel ca ale Pick a_uniform random variable to 
etermine the type failure. Decrement the 

appropriate stock. 


Call QUEUE: Enter queue with failed part. 


Check to see if server is available and the queue 
size is less than one: 


Is SERVR = 0 and QSIZE = 0 ? I£' so, enter service, 
go to step 15, else join queue, go to step 14. 


Part joins queue at the end of the service 
line. Increment QSIZE, store QPART, QTIME, check 
QMAX versus QSIZE, go to step 7. 


Enter service, store waiting time in queue _. 
(WAIT(i)), part type (QPART(1) ) advance positions 
ils the queue up one position. Galculate service 

ime: 

SERVTM = - ln (U) / SRATE(i) 
Calculate completion of service time (COMPLT): 
COMPLT = MSNTIM + SERVTM 


Store the stock levels at the moment of a given 
failure (JFAIL) ata given time (MSNTIM). Store in 
LOWLVL the lowest stock level at that time. 


If stocks are depleted, cease failures until a 
Eebadr OF the proper type occurs. FLAG2 = 0, Go to 
step : 


Call ALTQ1: Talley waiting time, queue size of. 
those parts which remain in the queue at mission 
termination time MSNTIM 


Is this run the last run: RUN s¢ TOTRUN ? 
If not, go to step 5 


Call STATS: compute statistics; 
Calculate mean up aircraft: MUPPLN = UPPLANE/RUN 


Calculate the variance and standard deviation of 
MUPPLN 


Calculate distribution of up aircraft; talley the 
number of up aircraft less than one standard | 
deviation from the mean, two standard deviations 
from the mean, or at the mean, etc CNTR). Talley 
the number of runs with zero aircraft up, one 
aircraft up, and so forth (YCNTR). 
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Step 22: Calculate the average wait in line WAITT, average 
maximum qsize MQBAR, and average number of 
customers waiting in the queue. 

Step 23: Print Statistics 

Dynamic 1 Model Modifications?*; 


ine. Are there equal numbers of parts ? I£ so, go 


Step 14A: een the number of parts of Oe i in_the queue 
to step 14c, else go to step 14b. 


Step 14b: Sort components by count, reorder position in line 
in favor of the most numerous component (to the 
front of the line). Return to step 7 


Step 14c: Sort parts by arrival time in the line (FIFO). 
Return to step 7. 


Dynamic 2 Model Modifications: 


Step 14a: Count the number of parts of type iin the queue 
line (COUNT(i)). 


Step 14b: Assign a priority to each type component i: 
RHO(i) = LAMBDA(i) / SRATE(i) 
LINEUP(i) = RHO(i) x COUNT(i) 

Step 14c: Sort components Sd lad Ae Eon 


line in favor of largest LINEU 
step 7 


er position in 
(i). Return to 
Dynamic 3 Model Modifications: 


Step 14a: Count the number of operating parts of each 
component i (STOCK(i)). 


Step 14b: Assign a priority to each component type based on 
the count of current stock in step l4a: 


LINEUP(i) < LINEUP(j) iff 
STOCK(i) < STOCK(3) 
Step 14c: Sort the parts in line by LINEUP(i reorder 
a P ~ t LENE ) 


poset oie in favor of the smalles UP (i 
eturn to step 7. 


FRL-L(H) Model Modifications: 


Step 2b: Call LINEUP; Read _in the values of LAMBDA(i), 
assign a priority from one to_n total type 
components based on the smallest (largest) 
LAMBDA(i). Return to step 3 


Step 14a: Sort que line by priority values assigned in 
step Zb. Reorder position in line to favor the 
eae one vcs with the smallest LAMBDA(i). Return to 
step 7. 


+*Modifications apply to FIFO algorithm, steps follow 
sequentially from there 
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APPENDIX D 
tab le ei. BRO URormiran Coding 


INITIALIZATION OF VALUES 


INTEGER RUN, TOTRUN, TOTSPR,FPART,QMAXMN( 500), 
TENTIV CNT 2, Cites, CNiteereNi dt, VENT 2, YCNTS,YCNT9, YCONT5,YCNI6, YCNT7, 
PVeENtopercino, VONILOPYeENT LI, YCNTI2, VENTS 


REAL PLANES, LAMBDAC10),SPAREC10),STOCK(C10),UPLAN( 500) 
1AC100000),LOWLVL(500),TOTALC10,500),AAC100000),CWAITC 
1PROBC10),MSNTIM, MQBAR, QT(500,120),FACTORC10), 
L1PPLANE,AAAC60000),MUPLN,SRATEC10) 


,LQ,K, 
100000), 


COMMON/ QUEU/STOCK,MSNTIM,L RUN, QMAXMN, LINEK, CWAIT,CUSTMR, REPARE, 
LISYSARR,TZZ,QT,SRATE,TOTMSN, COMPLT,“STAT/MUPLN,VRPLN,SDPLN,CNT1, 
1CNT2,CNT3,CNT4, YONT1, YONT2,YCNT3, YCNT4,YCNT5S,YCNT6,YCNT7, 
LYCNT8,YCNT9,YCNT10,YCNT11,YCNT12,YCNT13, 

IMQBAR,WAITT,LQ 


DATA JFAIL/0/,1/0/,11/07,RUN/1/7, TOTRUN/500/7, 
MUOUER U7 O7 pF EAGZZLZ, K/=5 .57 
TOTMSN = 120 
TOTSPR = O1 
PLANES =06.0 
READ IN FAILURE AND REPAIR RATES 

DORI = 2, TOLER 

READ (2,2) LAMBDACJ) 

FORMATCF8.5) 
CONTINUE 
DODS7) = 1, TOrPkT 

READ (3,8) SRATECJ) 

FORMATCF8.5) 
CONTINUE 
SUMLAM = 0 
DO 7 J = 1,7OTPRT 

SUMLAM = SUMLAM + LAMBDACJ) 
CONTINUE 

DO 9 I = 1,TOTPRT 
FACTORCI) = PLANES ® LAMBDACI) ® TOTMSN 
SPARECI) = NINTCFACTORCI)D + K ® SQRTCFACTORCI))) 
IFCSPARECI) .LT. 0) SPARECI) = 0 

CONTINUE 


CALL LLRAND II, GENERATE RANDOM NUMBER ARRAYS 
CALL RANDOMCA, AA, AAA) 
INITIALIZE A NEW MISSION 


JFAIL = 0 
PPLANE = PLANES 
XXX = SUMLAM * PLANES 
LRUN = RUN 
FLAG2 iN 
FLAG3 1 
DO 6 J = 1,TOTPRT 
STOCKCJ) = SPARECJ) + PLANES 
CONTINUE 


oat 


pel 


10 


20 


~> 


26 


BR) K MK 


21 
29 


30 


MK OK WG 
p— 


REPARE = 0 


CONTINUE 
IFCFLAG2 {CEQ 1) iien 


GENERATE ARRIVALS OF FAILED PARTS 


l= ie 
Z= I 
U = ACT) 
T = -LOGC U ) 7 © SUMLAM * PLANES) 
FAILT = T + MSNTIM 
ENDIF 
TEST TO SEE IF A FAILURE HAS OCCURED 
CONTINUE 
IFCREPARE .GE- 1] .AND. FEAG2Z “C771? THen 
REPARE = 0 
DO 25.452 ots Oeil 
IFCSTOCK(J) .LT. 1) THEN 
FLAG2 = 0 
FAILT = TOTMSN 
GO TO 2 
EES 
FLAG2 = ] 
ENDIF 
CONTINUE 
GO TO 10 
ENDIF 
REPARE = 0 


DO 26 J = 1,TOTPRT 
IFCSTOCK(J) .LT. 1) THEN 
2 


FLAG2 = 0 
FAILT = TOTMSN 
GO TO 23 
ELS 
FLAG2 = ] 
ENDIF 
CONTINUE 


INCREMENT TIME STEP 
VERIFY MISSION LENGTH 
IFCFAILLET LT. _COMPLT [AND SEE AGZ SCE aero ithen 
MSNTIM = FAILT 
IFCMSNTIM .GE. TOTMSN) GO TO 22 
GO TO 30 
EESE 
MSNTIM = COMPLT 
IFCMSNTIM .GE. TOTMSN) GO TO 22 
GO TO 2l 
ENDIF 


CALL ALTQ 
GO TO 20 
CALCULATE WHICH COMPONENT FAILED 


CONTINUE 


DO 31 J = 1,TOTPRT 
PROBCJ) = LAMBDACJ) 7 SUMLAM 
CONTINUE 
CALL A RANDOM NUMBER, USE LAMBDA'S TO BUILD CDF, 
AND DETERMINE WHICH PART FAILED 


_ 
ri 
ou 
1 
_—ri 
ab 
p= 


awn 
bh 


a2 
33 


34 
40 


4] 
G2 
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xX = AACTI) 
XLHR = 0.0 
VO ec ay Pe eay ci5 
XUPPR = XLHR + 
Te CCG eM IReeeaAND. X ~LE. AURPPR) THEN 
STOCNCI =) o OUCKCI) — 1 


FPART = J 
GO TO 33 
E 
XLWR = XUPPR 
ENDIF 
CONTINUE 


CALL QUEUECFPART, AAA, ¥34,%29,%70) 


CHECK STOCK LEVELS AND BRANCH 


CONTINUE 
ARE STOCKS LEVELS = 0? 
CONTINUE 
DO 41 J = 1,TOTPRT 
PeCSTOCKCJ): .GE. I) THEN 
FLAG2 l 


0 
TOTMSN 
Z 


G) 
© 
— 
© 
Hu il 


ENDIF 
CONTINUE 


JFAIL = JFAIL + 1 
RECORD STOCK LEVELS FOR A GIVEN FAILURE 


DO 51 J = 1,TOTPRT 
TOTALCJ,JFAIL) = STOCK(J) 
CONTINUE 


GO TO 10 


MISSION TERMINATION 

CALL ALTQI 
CONTINUE 

DOP ie = A 

LOWLVL (JJ 

L 


J).GT.TOTALCLL,JJ)) LOWLVLOJJ) = TOTALCLL,JJ) 
CONTINUE 
CONTINUE 


IFCLOWLVLCJFAIL) .LE. PLANES ) THEN 
UPLANCRUN) = LOWLVLCJFAIL) 

EESE 
UPLANCRUN) = PLANES 

ENDIF 


PLOT AND DISPLAY OUTPUT 


IFCRUN .GE. TOTRUND THEN 
WRITECO6,%) "CONTINUOUS TIME, CONSTANT FAILURE RATE, FIFO! 
fine C06,%) “MISSION TIME =",MSNTIMN, “JEAIL*,JFAILL 
WRITECO6,*) "INITIAL PLANES =',PLANES,* INITIAL SPARE FACTOR',K 
WRITECO6,*) "RUN NUMBER ',RUN,*CUMULATIVE FAIURES (CI)=',22 


(Nn 
LJ 


mx «| 


103 


1o3 


Z, £2 , TZ2 TRACK ARRAY INPUTS TO VERIFY THEIR DIMENSIONS ARE 
NOT EXCEEDED 


WRITECO6,%) ‘TOT PARTS=",TOTPRI,‘'2=" 34, 222 
HMRITECO6,%) 'LAMBDA=',XXX, 'SERVICE RATE= VARIABLE' 
WRITEC06,*) ‘TRAFFIC INTENSITY RHO= VARIABLE! 
ENDIF 
TEST RUN NUMBER TO STOP SIMULATION 


IFCRUN .GE. TOTRUN) GO TO 200 
RUN = RUN + 1] 
GO 1G 
COMPUTE SAMPLE STATISTICS 


CONTINUE 
CALL STATSCCUSTMR,LRUN, UPLAN, QMAXMN, LINEK, CWAIT, 
1TOTMSN, SYSARR, QT) 

WRITE OUTPUT OF STATISTICS 


WRITEC06,102) MUPLN,VRPLN,SDPLN 
FORMATC1X, ‘EXPECTED NUMBER OF UP PLANES =!,F11.8, 
17,1X,'VARIANCE OF UP PLANES =!,F15.4, 


17,1X,'STD. DEV. UP PLANES =", F15.9) 

WRITEC06,103) CNT1, 

1CNT2,CNT3,CNT4,RUN, TOTMSN, 

1,YCNT1,YCNT2, YCNT3,YCNT9¢,YCNT5,YCNT6,YCNT7, 
LYCNT8&,YCNT9,YCNT10,YCNT11,YCNT12, YCNT13 

FORMATC1X, *NUMBER OF RUNS UP PLANES <= -1 SIGMA =',14, 
17,1X, ‘NUMBER OF RUNS UP PLANES <= MEAN ='!,1G, 

17,1X, ‘NUMBER OF RUNS UP PLANES <= +1 SIGMA ='!,14, 
17,1X, ‘NUMBER OF RUNS UP PLANES <= +2 SIGMA =!,14, 
17,1X, "TOTAL RUNS =",15,"° “MISSION TInNeG=" hod, 


17,1X,*NUMBER OF RUNS W/Z O PLANES UP =!,15, 
17,1X,*'NUMBER OF RUNS W7 1 PLANE UP =!,15, 
17,1X, *NUMBER OF RUNS W7 2 PLANE UP =!,15, 
17,1X, "NUMBER OF RUNS bi7 3 PLANE UP =!,15, 
17,1X,*NUMBER OF RUNS W7 4 PLANE UP ='!,15, 
17,1X,*NUMBER OF RUNS W7 5 PLANE UP =',15, 
17,1X, ‘NUMBER OF RUNS W7 6 PLANE UP ='!,15, 
17,1X, "NUMBER OF RUNS W/ 7 PLANE UP ='!,15, 
17,1X, ‘NUMBER OF RUNS W7 8 PLANE UP ='!,15, 
17,1X, "NUMBER OF RUNS W7 9 PLANE UP =',15, 
17,1X,*NUMBER OF RUNS W410 PLANE UP =',15, 
17,1X, "NUMBER OF RUNS W711 PLANE UP =!,15, 
17,1X,'NUMBER OF RUNS W412 PLANE UP =!,15) 


WRITEC06,104¢) MQBAR,WAITT,LQ 

FORMATC1X,'AVE MAX QUEUE SIZE =',F5.0, 

17,1X,*AVE WAIT IN THE QUEUE =',F10.4, 

17,1X,*AVE NUMBER OF CUSTOMERS WAITING IN Q=',F10.9) 


DO 199 I = 1, 70OTPRI 

PRINT *,*SPARES TYPEC',1I,')=',SPARECI) 
CONTINUE 

PRIM, K-"5K 


STOP 
END 


Re CAA KA eK AK KK KK KK KR KK KO KARA RR RK KOH KM RRR RK MRK UNH 


SUBROUTINE QUEUECFPT,QA,x*,xX, *) 


SAVE QSIZE,QPART,LINEUP,SERVR,SERVPT,SERVIM,COMPLT, QMAX,RUN, II, 
LARRIVE,WAIT,QTIME 


INTEGER QSIZE,SPART, QPARTC200),SERVPT, 
1QMAX€500),QRUN,RUN,CNTR,FPT 


REAL QAC60000),QSTOCK(10),QT(500,120),WAIT(100000),QTIME( 200), 
ISRATEC10) 


COMMON /QUEU/QSTOCK, QMSNTM, QRUN, QMAX,CNTR,WAIT,QCUST,REPAIR, 
ISYSARR,2Z2ZZ,QT,SRATE, TOTMSN, COMPLT 


DATA SERVR/07,QSIZE/07,SERVPT707,14J707,RUN/0/4, 11707, 
IFLAG/0/ 


2 222. 222.20, 5,00, ,0,0,0, 2,0, 9, 0,8,9, 5, 4.0.9, 5,8. 5, 4,8. 4.8.8 tatatatetatetetatetetatatatatatatatetatatetatatetetetetatatatetetatetetates 


R 


PARTS ARRIVE 


TRGURUN Vii. 2 .AND. FLAG (tie 1) THEN 
FLAG = l 
QCUST = 0 
SYSARR 0 
DEPART 
COnP Ey 
ENDIF 


CLOCK = QMSNTM 
ARRIVE = 1 
SYSARR== SYSARR <t 1 


INITIALIZE QUEUE FOR A NEW RUN 
IFCQRUN .GT. RUN) THEN 
QMAXCQRUN) = 0 


0 
TOTMSN 


QSIZE = 0 

RUN = QRUN 

SERVR = 0 

COMPLT = TOTMSN 
ENDIF 


TS PHESSERVERS BUST? 
MECoeRVR «LT. TeeAND. QSI1ZE alu. 1) THEN 
QPARTC1) = FPT 
GO TO 10 
ENDIF 
GO TO 30 


CHECK SERVER AVAILABLE, QUEUE, AND SERVICE TIME 


ENTRY ALTQ 
ARRIVE = OG 
CLOCK = QMSNTM 


INITIALIZE QUEUE FOR A NEW RUN 
IFCQRUN .GT. RUN) THEN 
QMAXCQRUN) = 0 


WSIZE = 0 

RUN = QRUN 

SERVR = 0 

COMPLT = TOTMSN 
EW DIF 
Me@soeRVR .LT:.1 «AND, QSIZE GEEZ a1]) GO TO 10 
GO TO 100 


in 
cn 


30 


40 


10 


50 


100 


THE SQUiEVIES 


QCUST = CCuUsST + = 

IFCQSIZE =. 1) 3G0 Toma o 

QS TZE = OSes ae 

QPARTCQSIZE) = FPT 

QTIMECQSIZE) = CLOCK 

rena ALT. QSIZE) QMAXCQRUN) = QSIZE 


CONTINUE 

QSIZE = | 

QPARTC1) = FPT 

QTIMECQSIZE) = CLOCK 

toma -LT. QSIZE) QMAXCQRUN)D= 1 
0 


THE SERVICE AREA 


CONTINUE 
SPART = QPARTC1) 
STORE WAITING TIMES FOR QUEUE 


IFCQSIZE Gi. 0) THEN 

= eee 

WAITCII) = CLOCK - QTIMEC1) 
Sonne STORE NUMBER OF PEOPLE IN LINE IN VARIABLE CNTR 
N 


CNTR = II 


IF CQSTZE GT 221) hen 

PAIRS == eGolZe =o] 

DO 50 I = 1,PAIRS 
QPARTCI) = QPARTCI+1) 
QTIMECI) = QTIMECI+1) 


CONTINUE 
QPARTCQSIZE) = 0 
QTIMECQSIZE) = 0 
ENDIF 


IFPCQSIZE {Gl 2 O08 G51ZE = Gaze! 


SERVR = 1 
SERVE) = "SERVET er 1 


SERVICE IS EXPONENTIAL 


IJ =IJ+1 
XX = QACTIJ) 
Z22Z = IJ 


SERVTM = -LOGOCXX) 7 SRATECSPART) 
COMPLT = SERVI + Choe 
CHECK TO SEE IF SERVICE HAs “CGMP eeleD 
IF SO INCREMENT PART 


TLECCOMPET SEEs CEOCK)® THEN 


COMPLT = TOTMSN 
REPAIR = 1 
SERVR = 0 


DEPART = DEPART + 1 
QSTOCKCSPART) = QSTOCKCSPART) + 1 
ENDIF 


x RETURN TO APPROPRIATE SEGMENT OF MAIN PROGRAM 
re PEGREE Re EQ) LO CAND] QSIZE .GT.0) THEN 
REPAIR = 2 
GO TO 10 
CNDIF 


QTCRUN,CLOCK) = QSIZE 
Pe CAnine ve. Gl wl) ORETURN: 1 
RETURN 2 


x COUNT RESIDUAL PARTS IN QUEUE AT MISSION'S END 
ENTRY ALTQ] 
DEevolZe ie Gl. 0) (THEN 
DOW 29h = 1 @SIZE 


Pie Tite] 
CNTR = ITI 
WAITCII)D = TOTMSN - QTIMECI) 
129 CONTINUE 
ENDIF 
RETURN 3 
END 


KKK KKK KK KKK KKK KK KK KK KKK HK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KK KKK KK KK KKK KKKKRR 


SUBROUTINE STATSCQCUST,RRUN,UPPLAN, QMAX,LINE,SWAIT, 
LTOTMSN,SYSARR, QT) 


COMMON/STAT/MUPPLN, VARPLN,SDPLAN,CNTR1,CNTR2, 

ICNTR3,CNTRG, YONTRI,YCNTR2Z, YCNTR3,YCNTRG,YCNTRS, YCNTR6, YCNTR7, 
LYCNTR8, YCNTRO,YCNT10,YCNT11,YCNT12,YCNT13, 

IMXQBAR,WAIT,LQ 


INTEGER RRUN,QMAX(500), 
ICNTRI,CNTR2,CNTR3,CNTRG, YONTRI, YCNTR2, YCNTR3, YCNTRG, 
LYCNTR5,YCNTR6,YCNTR7, 
LYCNTR8,YCNTRI,YCNT10,YCNT11,YCNT12, YCNT13 


REAL MUPPLN,PSUM(C500),UPPLAN(C500) ,MXQBAR,SWAIT(100000), 


LTOTMSN,LQ,QTC500,120) 
HH HE HHH HE HEHEHE HH HEHEHE HE HEHEHE HH HH HEH HE HEHE HE HEHE HE HEH HE HE HEH HE HEH HK HH HEH HH HEH HH HH HEH HHH HH HK KKK KKK KKK 


WAIT = 0 
CNTR1 
CNTR2 
CNTR3 
CNTR 
YCNTR1 


YCNTR2 
YCNTR3 
YCNTRSG 
YCNTRS5 
YCNTR6 
TCHiR 7 
YCNTR8 
YCNTR9 
YCNT10= 0 
YCNT11= Q 
0 
0 


oooe 


ob 0d at heat 
coooo0eococo © 


YCNT12= 
YCNT13= 
MUPPLN 
VARPLN 
SDPLAN 
MXQBAR 
LQ = 0 
UPLANE 


Woy) tt 
© oooeo 


57 


105 


108 


MEAN CALCULATIONS 


SUM UP A/C AT MSN TERMINATION FOR ALL RUNS 
SUM UP QUEUE MAX'S FOR EACH RUN 


DO 104 J = 1,RRUN 
UPLANE = UPLANE + UPPLANCJ) 
Q = Q@ + QNAXCJ) 

CONTINUE 

IFCUPLANE =Ci00 THEN 
MUPPLN = UPLANE/RRUN 


EES 
MUPPLN = 0 
ENDIE 
DO 105 J = 1,RRUN 
PSUMCJ) = CMUPPLN - UPPLANCJ) ) x2 


SSPLAN = SSPLAN + PSUMCJ) 
CONTINUE 


MEAN MAX Q SIZE FOR ALL RUNS 
MXQBAR = Q 7 RRUN 


PRINT *,*LINE',LINE, *MXQBAR*,MXQBAR, 'Q',Q 
CALCULATE THE AVE. WAIT IN QUEUE 


DO 108 IL = 1,LINE 

WAIT = WAIT + SWAITCIL) 
CONTINUE 
WAIT = WAIT 7 SYSARR 


ASSOCIATED VARIANCE/STD DEV CALC'S 


IFCSSPLAN .GT. 0 .AND. RRUN .GT.1) THEN 
VARPLN = SSPLAN 7 CRRUN - 1) 
SDPLAN = VARPLNX*.5 


Ord 
0.0 


SDPLAN 


= MUPPLN + SDPLAN 

= MUPPLN + C2*SDPLAN) 

106 J = 1,RRUN 

se COC antou sLE. Xl) CNUIRI =" CNTR + 1 
1 

Jf 


~< 
_— 
| 
= 
vU 
“U 
— 
= 
1 


FCUPPLANCJ).LE.X2 .AND. UPPLANCJ).GT.X1) CNTR2 
FCUPPLANCJ).LE.X3 .AND. UPPLANCJ).GT.X2) CNITR3 
FCUPPLANCJ).LE.X@ .AND. UPPLAN(CJ).GT.X3) CNIRG 


CONTINUE 


Ble 


CNTR2 + 1 
CNTR3 + 1 
CNTRG + 1 


DO 107 J = 1,RRUN 


PEGUgeeeit ene.) Od YCHTRI = YCNTRI + 1 
PeCUmee Joe bee SAND. UPPLANCJ).GI.0) YCNTR2 = YCNIR2 + 1 
CUP een eee 2a. UPPLANCJ).GI.-LRYCNTIRS = YCHTRS + 1 
Petre tee Sees UPPLANCJ) -G).2) YOMERG = YCHTIRG + 1 
PeCUr Bets) Le eae AND UPPLAHCJ) °GT.3) YCMIRS = YCNTRS + 1] 
Pe CUR; en co eee AND. UPPUEATCJIOGT.49) YOMTR6 = YCHTR6 + |] 
PeOUmr anes eno ANDESUPPLAHCJJ°GI.5) YCHTR?7 = YCNTR7 + I 
PeaGueweanC ote. ? 6. AND. UPPLANCJ).GT.6) YCNTR&S = YENIR& + 1 
TeGUBp onic J be. 6 GAND. UPPLANCJ).GI.7) YCHTR9 = YENTR9 + 1 
TCU peewee ee. 9 SAND UPPLANCJ).GI.8) YCNTIO= YCNTIO+ 1 
TROUPE iCweee LOC AMD. UPPEAHCJ).GI.9) YCHTITI= YENTIi+ 1 
IFCUPPLANCJ).LE.11.AND. UPPLANCJ).GT.10) YCNT12= YCNT12+ 1 
PeOUrr CAG) LE ste AND URrEANs soe ol. Ll) YENTIS= YCNT13+ 1 
107 CONTINUE 
TEMP = 0 
DO 208 I = 1,RRUN 
DO 209 J =1,TOTMSN 
TEMP s= TEhin ste acl, J) 
209 CONTINUE 
208 CONTINUE 
LQ = TEMP 7 CRRUN X* TOTMSN) 
RETURN 
END 


HEHE HEHE I HH HEHEHE HE HE HE HE HEHE HE HE HE HE HE HEH HE HE HE HE HE HE I HE HE HE HE HE HE HE HE HE HE HE IE HE HE HE I HE HE HEH HE HE HK HEH HH HE HEH HH HH HE | 
SUBROUTINE RANDOMCRA, RAA, RAAA) 
REAL RAC100000),RAAC100000),RAAAC60000) 


N = 100000 
NN = 60000 
ke ey Mae rest, 
IIX = 14999 
IlTIxX = 13999 


CALL LRNDCIX,RA,N,1,90) 

CALL LRNDCIIX,RAA,N,1,90) 
CALL LRNDCIIIX,RAAA,NN,1,90) 
RETURN 

END 


50 


Sal 


Sifagl 
32 


341 


37 


S32 
S51 


55 


Dynamic Modifications 


THE QUEUES 


QCUST + 1 
IFCQSIZE -Upee 2 eG ieee 
STORE DATA ON NEWLY 
QSIZE. = QSUZEe ar 
QPARTCQSIZE) = FPT 
QTIMECQSIZE) = CLOCK 
IFCQMAXCQRUN)D LT. 


QCUST = 


QSIZE) QMAXCQRUN) = 


ARRIVED PART 


QSIZE 


ASSIGN Q@ PRIORITY BASED ON NUMBERS IN THE Q 
INITIALIZE AND COUNT NUMBERS OF PARTS IN Q 


DO 31 I = 1,TOTPRT 
QCNTR(I)D = 0 
CONTINUE 


DO 32 1 = 2, 01k 
DO 321 J 7] esIZe 
IFCI .EQ. QPARTCJ)) 
QCNTRCI) = 
ENDIF 
CONTINUE 
CONTINUE 


DO 341 I = 1,TOTPRT-1 
IFCQCNTRCI) 
FFLAG = 9 
EESE 
FFLAG = 1 
60 10-37 
ENDIF 
CONTINUE 


THEN 


QCNTRCI)D + 1 


~EQ. QCNTRCI+1)) THEN 


SORT AND ASSIGN PRIORTIES OF SERVICE 


IFCERELAG .<GE. JOS nREN 


DO 331 I = 1,TOTPRI 
DO 332 J = 1,QSIZE 
IFCQPART(J) .EQ. I) THEN 
LINEUP(J) = QCNTRCI) 
ENDIF 
CONTINUE 
CONTINUE 


ASSIGN PRIORITY TO THE LINEUP WAITING IN THE Q 
SORT AND ASSIGN ACTUAL POSITIONS IN Q BASED ON PRIORITY 


PAIRS = QSIZE - 1 
DONE = 1 
IFCDONE .EQ. 1) THEN 
DONE = 0 
DO 351 I = 1,PAIRS 
IFCLINEUPCI)D . 
TEMPL1 
hEneee 
TEMPUS 
LINEUPCTI) 
QPARTCI 
QTIMECI 
LINEUP ( 
QPARTCI 
QTIMECT 
DONE = l 
ENDIF 


I 
) 
) 
i+ 
+1 
+] 


LT. LINEUPCI+1)) THEN 
LINEUPCI) 

QPARTCI) 

QTIMECI) 

= LINEUPCI+1) 

= QPARTCI+1) 

= QTIMECI+41) 

1) = TEMPL1 

) TEMPL2 

) TEMPLS 


S51 


36 


361 


40 


30 


oy 


CONTINUE 
PAIRS = PAIRS - 1 
COF 1035 
eB 0BUs 
ENDER 


TRECERCAG ee CO THEN 
PeniRow= VSIZE = 1 
PDONE = 1 
Peer OOKE. EQ.) ) THEN 
PDONE = 0 
DO 361 I = 1,PPAIRS 
PecQuluecr)” .Gl. QTIMECI+1)) THEN 


PTEMP2 = QPARTCTI) 
PTEMP3 = QTIMECT) 
QPARTCI) = QPARTCI+1) 
QTIMECI)D = QTIMECI+1) 
QPARTCI+1) = PTEMP2 
QTIMECI +1). = PIEMPS 
PDONE = 1 
ENDIF 
CONTINUE 
PPAIRS = PPAIRS - 1 
GO TO 36 
ENDIF 
ENDIF 
GO TO l 
CONTINUE 
QSIZE = 1 


QPARTC1) = FPT 
QTIMECQSIZE) = CLOCK 
Bae -LT. QSIZE) QMAXCQRUND= 1 


Dynamic 2 Modifications 


THE QUEUES 


QCUST = QCUST + 1 
IFCQSIZE .LT. 1) GO TO 40 
STORE DATA ON NEWLY ARRIVED PART 
QSIZE = QSIZE + 1 
QPARTCQSIZE) = FPT 
QTIMECQSIZE) = CLOCK 
IFCQMAXCQRUN) .LT. QSIZE) QMAXCQRUN)D = QSIZE 


ASSIGN Q PRIORITY BASED ON NUMBERS IN THE Q 


INITIALIZE AND COUNT NUMBERS OF PARTS IN Q 
DO 31 I = 1,TOTPRT 
QCNTRCI) = 0 
CONTINUE 


DO 32 I = 1,TOTPRT 
DOzSs21 J = 1,Q51ZE 
IFCI .EQ. QPARTCJ)) THEN 
QCNTRCI)D = QCNTROCI)D + 1 
ENDIF 
CONTINUE 
CONTINUE 


x SORT AND ASSIGN PRIORTIES OF SERVICE 


DO 331 I = 1,TOTPRT 
DO 332 J = 1,QSIZE 
IFCQPARTCJ) .EQ. I) THEN 
LINEUPCJ) = QCNTRCI)*RHOCI) 


ENDIF 

ge CONTINUE 
po) CONTINUE 
x ASSIGN PRIORITY TO THE LINEUP WAITING IN THE Q 
x SORT AND ASSIGN ACTUAL POSITIONS IN Q BASED ON PRIORITY 

PAIRS =3051ZE5 se 

DONE = ] 
aD IFCDONE .EQ. 1) THEN 

DONE = 0 


DO 351 I = 1,PAIRS 
IFCLINEUPCI) .LT. LINEUPCI+1)) THEN 


TEMPL1 = LINEUPCT) 
TEMPL2 = QPARTCTI) 
TEMPL3 = QTIMECI) 
LINEUPCI) = LINEUPCI+1) 
QPARTCI) = QPARTCI+1) 
QTIMECI) = QTIMECI+1) 
LINEUPCI+1) = TEMPLI 
QPARTCI+1) = TEMPL2 
QTIMECI+1) = TEMPLS 
DONE = 1 
ENDIF 
oa CONTINUE 
PAIRS = PAIRS =72 
GO TO 35 
ENDIF 
GO 10 1 
40 CONTINUE 
QSIZE = 1 


QPARTC1) = FPT 

QTIMECQSIZE) = CLOCK 

IFCQMAXCQRUN) .LT. QSIZE) QMAXCQRUND= 1 
GO TO l 


Cc. Dynamic 3 Modifications 


x THE QUEUES 


30 QCUST = QCUST + 1 
IFCQSIZE .LT. 1) GO TO 40 
x STORE DATA ON NEWLY ARRIVED PART 
QSIZE = QSIZE + 1 
QPARTC(QSIZE) = FPT 
QTIMECQSIZE) = CLOCK 
IFCQMAXCQRUN) .LT. QSIZE) QMAXCQRUN) = QSIZE 


62 


341 


SY) 


aac 
Ss 


35 


551 


36 


561 


COUNT NUMBERS OF OPERATING PARTS IN STOCK 


NOP al = 1 TOTERI >) 
IFCQSTOCKCI) .EQ. QSTOCKCI+1)) THEN 


FFLAG = 0 
ELSE 
FFLAG = 1 
GO! 7G 37 
ENDIF 
CONTINUE 


SORT AND ASSIGN PRIORTIES OF SERVICE 
TFCEELAG. GE. 4a) THER 


DO 331 I = 1,TOTPRT 
DO 332 J = 1,QSIZE 
IFCQPARTCJ) .EQ. I) THEN 
LINEUPCJ) = QSTOCKCI) 
ENDIF 
CONTINUE 
CONTINUE 


ASSIGN PRIORITY TO THE LINEUP WAITING IN THE Q 
SORT AND ASSIGN ACTUAL POSITIONS IN Q BASED ON PRIORITY 
PAIRS = QSIZE - 1 
DONE = 1 
1) THEN 


IFCDONE .EQ. 
DONE = 0 
DO 351 I = 1,PAIRS 
IFCLINEUPCI) .GT. LINEUPCI+1)) THEN 


TEMPL1 = LINEUPCI) 
TEMPL2 = QPARTCT) 
TEMPLES = QTIMECT) 
LINEUPCI) = LINEUPCI+1) 
QPARTCI) = QPARTCI+1) 
QOTIMECI) = QIIMECT SID 
LINEUPCI+1) = TEMPL1 
QPARTCI+1) = TEMPL2 
QTIMECI+1) = TEMPL3 
DONES= 1 
ENDIF 

CONTINUE 

PAIRS = PAIRS ~ 1 

GO TO 35 

ENDIF 


ENDIF 


PRCEEEAG «labo 0) THEN 
PPAIRS = QSIZE - 1 
PDONE = 1 
IFCPDONE .EQ. 1) THEN 
PDONE = 0 
DO 361 I = 1,PPAIRS 
IFCQTIMECI) .GT. QTIMECI+1)) THEN 


PTEMP2 = QPARTCT) 
PTEMP3 = QTIMECI) 
QPARTCI) = QPARTCI+1) 
QTIMECI) = QTIMECI+1) 
QPARTCI+1) = PTEMP2 
QTIMECI+1) = PTEMP3 
PDONE = 

ENDIF 


CONTINUE 


6 3 


40 


m— ROK 


WG 0° XK 


PPAIRS =PPATRS=—— 
GO TO 36 
ENDIF 
ENDIF 
GO TO 1 


CONTINUE 

QSIZE = 1 

QPARTC1) = FPT 

QTLIMECQSIZE) =—Chock 

coco a ee -LT. QSIZE) QMAX(CQRUN)= 1 


FRP-L(H) Modifications 


READ IN FAILURE RATES 
DO 1 J = 1,TOTPRT 
READ (2,2) LAMBDACS) 
LAMBDACJ) = 0.0167 
FORMATCF8.5) 
CONTINUE 
DO 3.3 = 1 TORR 
READ (3,8) SRATECJ) 
SRATEC JS) = 0.84 
FORMATCF8.5) 
CONTINUE 
SUMLAM = 0 
DO 7 J = 1,TOTPRT 
SUMLAM = SUMLAM + LAMBDACJ) 
CONTINUE 
DO 9 T=1,TOTPRT 
FACTOR(I) = PLANES * LAMBDA(I) %* TOTMSN 
SPAREC(I) = NINTCFACTOR(I) + K ® SQRTCFACTORCI))) 
IFCSPARECI) .LT. 0) SPARECI) = 
CONTINUE 


ASSIGN PRIORITY TO PARTS VIA SUBR. QPRIOR 
CALL LINEUPCLAMBDA,TOTPRT) 
CALL LLRAND II, GENERATE RANDOM NUMBER ARRAYS 


CALL RANDOMCA,AA,AAA) 
INITIALIZE A NEW MISSION 


2 2.2.2, 2,2,2,2,2,2,0,5,2,5,5,5,5,2,5, 5,0, 5,5,5,5, 5,45, 5,5,5,5,5,5, 5,5, 5, 5,5, 5,5,5, 5.5.5.5 05. Satetatetatetotetotatatetatotatatetetatetetatatete: 


% 


30 


505 
304 


THE QUEVES 


QCUST = QCUST + l 
IFCQSIZE .LT. 1) GO TO 40 
QSIZE = QSIZE + 1 
ele oa = ce 
IMEC QSIZ = 7 
TECQMAXCGRUN) JLT. QSIZE) QMAX(QRUN) = QSIZE 


DO 304 J = 1,TOTPRT 
DO 305 I = 1,QSIZE 
IFCQPART(I) .EQ. J) THEN 
LINEUP(I) = PRIOR(J) 
ENDIF 
CONTINUE 
CONTINUE 


301 


302 


505 


40 


aS = AS al 
DONE = 1 


A DOMWE Seo. Tarren 

DONE = 0 

DONS 0Zel = sl PAIRS 

PeCEINEUP GI). Gl2 JERR EUR Gil ale)) THEN 
TEipe= QTIMEGT) 
TEMPQ = QPARTCI) 
TEMPOYe="LIREUPC TD 
OI MECC =" OT IMECHE 
QPARTCI) = QPARTCI+1) 


BIREUPCID = LINEURGIA1) 
QTIMECI+1) = TEMP 
QPARTCI+1) = TEMPQ 
ETHEUPCI+1) = ,EMPOr 
DONE = 1 
ENDIF 
CONTINUE 
PAIRS = PAIRS - 1 
Go 10) 30) 
en DT 
DO 303 I = 1,QSIZE 
CONTINUE 
COenra 1 


CONTINUE 

QSIZE = 1 

QPARTC1) = FPT 

QTIMECQSIZE) = CLOCK 

2 aa -LT. QSIZE) QMAXCQRUN)D= 1 
G 1 


oo 2S Sete Pee eee see eee eee cee cece cece 0 ¢2 525.5 55255 55 5255.5. 5555.5, 5,5,5,5,5,2,5,5,5,25,5,5,5,3 


SUBROUTINE LINEUPCQLAMBD, QTOTPT) 
INTEGER QPRIOR(100) 

REAL QLAMBD(10),QTOTPT 

SAVE QPRIOR 

COMMON/RULE/QPRIOR 


DO 1 J=1,QTOTPT 
QPRIORGJ) = J 


CONTINUE 

PAIRS = QTOTPT - l 

DONE. = 1 

REGDONE .GE. 1JTREN 
DONE = 0 


DO 3 J=1,PAIRS 
IFCQLAMBD(J) .GT. QLAMBDCJ+1)) THEN 
TEMey =-QPRIGkt 
OPRIOR CUD -SOPk LOR GI2 .) 
QOPRIOKCS 21 = Genie 


DONE = l 
ENDIF 
CONTINUE 
PAIRS = PAIRS - 1 
GO" 10 2 
ENDIF 
RETURN 


END 
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